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TABLE OF SYMBOLS AND ABBREVIATIONS 


Amplitude of flow oscillation 
Fourier coefficient 

Added mass coefficient 

Drag coefficient 

RMS value of calculated force 
Inertia coefficient 

Cylinder diameter 

In-line force 

Measured in-line force 
Keulegan-Carpenter number, K = Ut/D = 20) A/D 
Critical Keulegan-Carpenter number 
Mean sand-roughness height 
Relative roughness 

Cylinder length 

Pressure 

Reynolds number, Re = U,D/v 
Period of flow oscillation 
Time 

Instantaneous velocity 
Maximum velocity in a cycle 
x-component of velocity 


V=component of velocity 


Horizontal component of coordinate axis 
Vertical component of coordinate axis 
Frequency parameter, 8 = Re/K = D?/vT 
Angle, 69=s2ne7 7 

Kinematic viscosity of water 

Density of water 


Circular frequency 
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I. INTRODUCTION 


The determination of the forces acting on bluff bodies 
immersed in time dependent flows has been of special 
interest in many branches of engineering. Special attention 
has Been devoted to the wave forces on offShore structures. 
The extensive laboratory and field studies of the past 
decade have yielded extremely interesting results regarding 
the in-line and transverse forces acting on smooth and rough 
cylinders [Ref. 1]. However useful these investigations 
were, they were based on the use of an empirical equation 
(Morison's) (Ref. 2] and on the force coefficients obtained 
either through the Fourier Analysis or the method of least 
Squares. In recent years, it became quite clear that 
further progress toward the understanding of the physics of 
the phenomenon must be based on analytical studies and the 
understanding of the kinematics of the flow. 

A theoretical analysis of the separated flow about a 
Circular cylinder is an extremely difficult problem. 
Efforts in that direction through the use of the discrete 
vortex models did not lead to meaningful conclusions. It 
became clear that the key towards the understanding of flow 
Separation, transition and vortex formation is in the very 


low Keulegan-Carpenter number (K) regime. It iS because of 
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this reason that experiments were carried out in an 
oscillating water tunnel for force measurementS and with an 
oscillating cylinder in a water table for flow visualiza- 
tion. The investigation resulted in findings beyond the 
original objectives. In the following, first the 


theoretical analysis and then the experimental apparatus and 


v 
¥ 


measurements will be presented. These will be followed by a 


discussion of the results and conclusions. 


sai 


II. ANALYSIS 


A. INTRODUCTION 

This section will present a condensed review of the 
theoretical studies regarding the oscillating flow about 
pie modes First the flow about an infinite flat plate 
undergoing sinusoidal oscillations will be briefly 
discussed. Then Schlichting's [{Refs. 3 and 4], Stokes' 
{[Ref. 5 and 6], and finally Wang's [{Ref. 7] solution of a 
sinusoidally oscillating cylinder in a still fluid will be 
discussed. It should be emphasized that the case of a 
cylinder oscillating in a fluid otherwise at rest is 
identical to that of a fluid oscillating about a cylinder. 
The only difference comes from the 'buoyant' force acting on 
the cylinder in the direction of flow due to the pressure 
gradient needed to accelerate the flow, i.e., the inertia 
coefficient for the case of the oscillating flow is equal to 
(1 + Cc.) where C. is the added mass coefficient for the 


oscillating cylinder, i.e., C = 1] + Ca: 


B. STOKES' OSCILLATING PLATE PROBLEM 

The flow about an infinite flat plate undergoing 
Sinusoidal oscillations was treated by Stokes assuming that 
the flow is laminar and stable. With respect to the 
coordinate axis in Figure 1, the Navier-Stokes equations and 
the boundary conditions reduce to: : 


Le 





du fe 
= > 5 (1) 
One dY 
ileal condi 10M: u(y,9) = @ 
boundary conditions: u(@,t) = U, cos wt 
u(~,t) = @ 


The solution of Eq. (1) via the Separation of variables 


is given by: 


u(y,t) = Uy exp(-y/>~ ) cos(ut - y/>~ ) (2) 


Eaeeing 7 = Via , one has 


u(y,t) = U. exp(-n) cos(wt - n) (3) 
"The velocity profile u(y,t) has the form of a damped 
harmonic oscillation, the amplitude of which 1s Ue exp (-y/=2 
in which a fluid layer at a distance y has a phase lag with 
respect to the motion of the wall [Ref. 3]." Figure 2 shows 

the velocity distribution in the neighborhood of the 


oscillating wall. 


eeeoCHLICHTING'S OSCILLATING CYLINDER PROBLEM 
Schlichting [Refs. 3 and 4] considered the case of a 


cylinder oscillating at high frequency and small amplitude 
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in a fluid otherwise at rest and used the boundary layer 


approximation to the Navier-Stokes equations. This lediicon 





2 
au au au _ _ 1 ap @“u 
at | ox ” “ ony e GocnlnaneeD a 
ay 
which reduces to: 
Bee 20. yg au 
a ax att YU 3x (3) 


for the outer flow, i.e., at the edge of the boundary layer. 


Combining Eqs. (4) and (5), one has: 





au ae 
1 ree A 0) aU 
It CO ee oe (6) 
ay 


This equation may be solved to the desired degree of 


approximation. The first order (linear) approximation 





yields: 
du o> 
1 Il _ au 
a a wm) 


For the approximation to be valid, the convective 
acceleration must be considered smaller than the local 


acceleration, i.e., 
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& 
@ 
Cc 

ras 

aN 
@ 
Cc 


Using the order of magnitude arguments, one has: 


3U dU 
» U ax” Dig Un” 


It is seen that for the analysis to be valid, one must have 


U_T 
m A . ; : : 
=~ << 1). O& D << 1 where Un is the maximum velocity in a 


cycle; T, the period of oscillation; D, the diameter of the 
cylinder; and A, the amplitude of oscillation. 

The solution of Eq. (7) with the following boundary 
Somailtions: 


aibary Detioms HED aco mare LO) 


it 
a] 
Cc 

i 


UN 
8 


at y , u= U 


yields the following stream function: 
LW 
v (Seyret) = f2 Uy(x) E(nper”™ (9) 
from which one obtains the velocity as: 


u(x,y,t) = U(x) [cos(ut) - exp(-n//2)cos(ut -nW/2)] 


(19) 


I 


which iS seen to be identical to Stokes' solution for the 
oscillating infinite flat plate. 
A second approximation may be applied to Eq. (4) by 


retaining the convective terms. Then one has: 





“= - SU Uo) (11) 


The stream function needed for the solution of the above 


equation is given by: 


du 
vy(eryety = JZ usc gtd (epp(me te + Eng(n)} (12) 


where £4, represents the periodic motion and f5,, the steady 


p 


State. The use of the following boundary conditions on Cop? 


at n= Qg u == @ ==> f= £' = Qg 


at On Sas Eo =26 
and on f5.: 


ac 0 = 26 lo =“7 =F0 => i eee 


at n= .o f' = @ 
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leads to an important result. The boundary condition at 
n=, £* = 6, for the steady state component cannot be met. 


Instead a finite value is obtained: 
me Be dau 
wes) = Soo 8 la (13) 


Periodically oscillating flow about a circular cylinder 
produces a steady streaming motion outside the boundary 
layer. In Figure 3 the streamlines about an oscillating 
cylinder ina still fluid are shown. This is for the case 
for a small Reynolds number with K -> @ which leads to four 
recirculation cells and four stagnatiton points around the 
cylinder. 

The force acting on the cylinder has not been evaluated 


by Schlichting. 


D. STOKES' OSCILLATING CYLINDER PROBLEM 

Stokes, as one of his considerable achievements, 
investigated the flow about oscillating cylinders and 
Spheres without resort to the boundary layer approximation 
(Refs. 5 and 6]. His solutions are then valid at an as yet 
unspecified small Reynolds number. His solution yields the 


following force per unit length of a cylinder: 








hy 


which can be obtained by use of the asymptotic expansion of 
the Bessel function for large values of 8, (8 = D*/vT) . 
Stokes' solution will be taken up again later following the 


discussion of amore recent solution. 


E. WANG'S OSCILLATING CYLINDER PROBLEM 

Wang [Ref. 7] in 1968, investigated the flow about an 
oscillating circular cylinder by taking into consideration 
the effects of curvature which were neglected in 
Schlichting's analysis. Using inner and outer expansion and 
the boundary layer approximation, Wang obtained the 
following force expression for the oscillating cylinder, 


valid to the ordermor g~3/2, 











2 2 
F 4n il it li T 4 l : 
ee ( + — - ———)coswut - -(1 + — + )sinwt 
=0U,,°D : / 18 ue 418/18 : [8B 1B / 1B 
+ higher order terms Cl 5} 


A comparison of Eqs. (14) and (15) shows that the first 
two terms in the corresponding parentheses are identical and 
Wang's solution differs from that of Stokes only in the 
third term. It is noted that both solutions are valid only 
for large 8 and thus they yield virtually identical results, 
Since the third term in each parenthesis in Eq. (15) is 


negligible. 
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F. MORISON'S EQUATION AND THE STOKES AND WANG SOLUTIONS 
The normalized form of the Morison equation, with 


U = U. cos wt, {Ref. 2] may be written as: 


— = Ca [cos 6] cos 6 - rg Sin @ (16) 
Beene 
which was developed for the determination of wave forces on 
offshore structures and discussed in the literature in great 
detail during the past 38 years [Ref. 1]. 
The force expression obtained by StokeS and Wang may be 
expressed in the form of Morison's equation by linearizing 


feos 6] cos 6. 


From (Ref. 8], writing: 


Bie 


ieas 6 | Cece ) g ek SuIme@os ce cos mmcma. (17) 
T 
n=@ 
ip cos? ng de 
—= == e©e°@ 
ay - ay cos 6 + a5 cos 286 
with: an = @ n even 
n+1 
A =i) é —_*—_ n odd 
‘ n(n” - 4)7 
eae: 
ai 3a 


Ie) 


using a £lrst term approxineeton. 
cos @/ cos 6 = 


cos 86 (533) 


Then Eq. (15) may be reduced to: 














i i 1 
C,=2+ 4 + ies (35) 
m xB 1B aR 
3n7 f 1 1 1 
Cy = x + i ad (20) 
[ 7B " 478 /1wB 
for the case of the oscillating flow about a cylinder 
(note: CH = l1 + C.)- These drag and inertia coefficients 


will be compared with those obtained experimentally with 


smooth and rough cylinders. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


Force meaSurementS were carried out in a large U-shaped 
water tunnel. It is 34 feet long, 25 feet high, and has a 
working cross section of 3 feet by 4.66 feet. Additional 
details of the design and operation of the water tunnel are 
given in [Refs. 9 and 18]. The flow in the tunnel is 
oscillated by means of a1 HP fan and a butterfly valve 
system. The frequency of oscillation of the butterfly valve 
is adjusted by means of an electronic feedback control 
system, so as to oscillate the flow at its natural 
frequency. The amplitude of oscillations can be maintained 
at the desired level by opening or closing a gate at the 
exit of the fan. 

The test cylinders, approximately 3 feet in length, are 
mounted horizontally and supported at both ends by force 
transducers. These transducers can measure both the in-line 
and transverse forces. The electronic signals from the 
force transducers and from the amplitude gages are fed to 
two amplifiers and then digitized by means of an A/D 
converter. The digitized signal is then analyzed by means 
of a computer. 

The experiments for each cylinder were repeated at least 


twice to ascertain the repeatability of the data. 


Px 


Two smooth cylinders with 8 = 18035 and 8 = 138@ and one 
rough cylinder with 8 = 1888 were tested. 

Flow visualization experiments were carried out ina 4 
foot by 8 foot water table. Circular cylinders of suitable 
diameter were oscillated by means of a slider-crank 
mechanism and a variable speed motor. The flow pattern 
about the cylinder was made visible by means of an electro- 
chemical technique [Ref. l1l]. For this purpose a solder 
wire was inbedded along the cylinder length. Voltage of 
5 - 18 volts were applied between the solder and an 
electrode (placed away from the cylinder) to generate a 
metallic white smoke from the solder. A small amount of 
NaCl was added to the water to increase conductivity. 
Careful adjustment of the applied voltage and the repeated 
cleaning of the solder and cylinder surfaces were necessary 
for the successful visualization of the resulting flow 


patterns. 
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IV. PRESENTATION AND DISCUSSION OF RESULTS 


The drag and inertia coefficients were evaluated through 


the use of Fourier averaging [Ref. 1], as given by: 


av 
: 3 en cos @ _ 
é | wcamapiaamaale 2 
; eDU, 
27 
OK ae Sin 980 
G = — de 
m 3 DD 
1 g o DU 


Figures 4 and 5 show the drag and inertia coefficient, for a 
Smooth cylinder, as a function of K for 8 = 10435 and B = 
13889. Similar data are shown in Figure 6 for a rough 
cylinder (K/D = 1/1988) for 8 = 18@@. Also shown in 

Figures 4 through 6 is the drag coefficient obtained 
Bueoretically by Wang, Eq. (15). 

Figures 4 and 5 show that Ca follows the theoretical 
line for a K value smaller than a critical Kar and then 
increases almost abruptly to a higher value. Subsequently, 
Ca decreases with increasing K and reaches a minimum at a K 
value of about 2. The variation of Cqg with K for K > 5 has 
been well documented in the literature [Ref. 3] and will not 


be dealt with here further. 
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Figure 5 is similar to Figure 4 in many respects with 
the addition of a hysteresis in the drag coefficient. For 
very small values of K, Ca followed the Wang-line and then 
rapidly jumped to a higher value. However, when the 
experiments were carried out by decreasing K at suitable 
intervals of K (with 15 minutes between each sampling), CE 
remained at its higher value and then rapidly jumped back to 
the Wang-line at a K value smaller than that found in 
experiments with increasing K. The reason for the observed 
hysteresis is not easy to explain and it may be due to the 
precarious nature of the stability of the flow to the 
prevailing disturbances. It is tempting to think that in 
the experiments with decreasing K the level of the ambient 
turbulence is higher and this forces K,, tO occur at a lower 
K. In the future, it may be desirable to wait longer 
periods of time between two successive changes of K. 

Figure 6 does not show any region where Ca agrees with 
Wang's prediction. In fact, from the lowest K obtainable in 
the present experiments to the K at which Ca goes through a 
minimum, Wang's prediction is considerably lower than that 
obtained experimentally. One will have to postulate that 
the flow has already become critical at K values smaller 
than about @.45 and one is observing only the post-critical 
region of the drag coefficient. 

Figures 4 through 6 show that the inertia coefficient 


obtained experimentally agrees quite well with that given 
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by Wang for values of K smaller than K,,- It is noted that 
both the theoretical and experimental values of Cn are 
Slightly larger than 2 due to the effects of Vso weve 1 
the boundary layer. It should also be noted that there is 
very little or no scatter in the C, data because the flow is 


in the inertia dominated regime. 


“y 


Figure 7 shows the RMS (root mean Square) value of the 
normalized in-line force as a function of K for all the 
cylinders tested. The theoretical value of the said 


coefficient is given by: 





which is seen to agree well with the data for K smaller than 
about 8. This is due to the fact that the effect of Cg on 
Ce (rms) is quite small and the value of C_e (rms) is dictated 
by Cae 

The in-line force measurements and the resulting drag 
and inertia coefficients have revealed certain heretofore 
unknown facts which are both interesting and surprising. 
Interesting in the sense that the flow exhibits numerous 
transitions as K varies from about 9.4 to say 290. 
Surprising in the sense that one would not have expected 


transitions and hysteresis at relatively small K values 
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particularly prior to the separation of flow which was 
expected to occur at a K value in the vicinity of 1.6 - 2.4. 
These two observations led to the flow visualization studies 
noted earlier. 

Figure 8 shows a sample photograph of the flow pattern 
about an oscillating cylinder. Clearly, Taylor-Gértler 
oe form around the cylinder at regular intervals when 
the amplitude of oscillation exceeds a particular value. 
Figure 9 shows a plot of K versus 8 for smooth cylinders. 
For K values smaller than a critical K and for a given 8B the 
flow remains stable and attached. The symbol 'O!' denotes 
the K value at which the Taylor-Gortler instability has 
manifested itself. This instability gradually leads to 
Separation and turbulence as K is increased. The symbol 'X' 
denotes the K values at which the flow has become turbulent 
in the boundary layer beyond which no clear vortices were 
observable. 

It 1s in light of Figure 9 that one can interpret the 
changes in Cq in Figures 4 through 6. For 8 = 1835 the 
instability starts at ees 1.85, according to Figure 9. In 
Figure 4 the jump in Cg occurs at a Slightly lower kK. | 
Nevertheless, it is now clear that the rapid changes in Cy 
in Figures 4 and 5 are a consequence of the instability 
observed in the flow. It is also seen that the K value at 


which the boundary layer becomes turbulent nearly 
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corresponds to the K value at which the drag coefficient 
goes through its minimum. 

The foregoing discussion has not addressed itself to the 
case of rough cylinders. Exploratory experiments with a 
sand roughened oscillating cylinder have shown that the 
effect of roughness is to reduce K,,- In fact, it was found 
that Ke = 9.4 at 8 = 1800. Such a low K value was not 
attalinabile in the force measurements with the rough 
cylinder due to a number of reasons, the most important 
being the noise level in the signal. Thus, it is safe to 
assume that the flow for the rough cylinder has become 
unstable at a K value at about 8.4 or less and led to the 
rapid increase in Cy prior to the lowest shown in Figure 6. 

In summary, it 1s now clear from the measurements and 
observations that the theoretical predictions of Stokes and 
Wang agree with the measurements as long as the flow is 
lmainar, attached, and stable. As the flow becomes unstable 
to Taylor-Gortler vortices, the drag coefficient is 
drastically affected primarily due to changes in the shear 
stress level and distribution over the cylinder. The next 
important change in the flow and in the force coefficients 
takes place when the flow separates and the boundary layer 
becomes turbulent. The inception of separation nearly 
corresponds to the occurrence of minimum drag. Beyond that 


the pressure drag becomes increasingly dominant and the 


Ze 


shear becomes negligible as K increases. Contrary to the 


previous assumptions, that Cg varies from 8 to Its maximum 
value aS K varies from 9 - 12, it is seen that there are 
both interesting and important phenomena in the low 
Keulegan-Carpenter number regime. The experiments described 
here, however difficult they may have been due to the 
problems associated with the measurements of very small 
forces, have shed considerable light on the low K flow 


regime. 
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V. CONCLUSIONS 


The investigation reported herein warranted the 


following conclusions: 


il. 


The theoretical solutions of the oscillating flow 
about a circular cylinder are valid only for attached, 
Stable, and laminar flow; 


The flow becomes unstable to Taylor-Gértler vortices 
at a critical Keulegan-Carpenter number for a given 
ratio of the Reynolds number to the Keulegan-Carpenter 
number, 1.@., 8; 


The drag coefficient shows a rapid jump at or near the 
Critical Keulegan-Carpenter number. This is often 
accompanied by hysteresis; 


There is a second critical Keulegan-Carpenter number 
at which the boundary layer separates and becomes 
turbulent. This paricular value of K nearly 
corresponds to the K value at which the minimum drag 
OCECULS ; 


The theoretical value of the inertia coefficient 
agrees with that determined experimentally only in the 
attached, laminar, and stable flow regime; and 


Roughness precipitates transition and increases the 
drag coefficient. The effect of roughness on the 
inertia coefficient in the stable flow regime is 
negligible. 
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Figure 1. Coordinate Axis for an Oscillating Flat Plate 
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Figure 2. Velocity Distribution for Stokes' Oscillating 
Flat Plate 
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